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I.  INTRODUCTION 


We  are  interested  in  determining  validated  sets  of  elementary  chemical 
reactions  for  use  in  predictive  combustion  models •  The  approach  used  has  been 
to  simulate  the  one  dimensional,  premixed,  laminar  steady-state  flame.  This 
approach  has  the  advantage  that  the  predicted  temperature  and  species 
profiles,  as  well  as  the  flame  speeds,  can  in  principle  be  compared  with 
suitable  burner  experiments  of  the  flame. 

Such  a  model  requires  as  input  data  sets  of  (1)  thermodynamic 
coefficients,  (2)  transport  coefficients,  and  (3)  chemical  reactions, 
including  their  rate  coefficients.  Often,  many  of  the  input  parameters,  and 
in  particular  the  rate  coefficients,  are  not  well  known.  Therefore,  it  is 
useful  to  be  able  to  perform  a  sensitivity  analysis,  so  as  to  systematically 
determine  the  effect  of  uncertain  parameters  on  the  solutions  of  the  model. 

If  the  system  is  very  sensitive  to  a  given  parameter,  more  time  and  effort  may 
be  justified  to  determine  a  more  accurate  value  for  that  parameter.  On  the 
other  hand,  an  insensitive  parameter  requires  less  effort;  and  it  may  be 
possible  to  eliminate  it  and  thus  simplify  a  complex  system.  In  addition,  the 
sensitivity  analysis  is  useful  in  understanding  a  complex  mechanism,  since  it 
indicates  which  parts  of  the  mechanism  are  important  for  a  given  problem. 

We  are  also  interested  in  the  region  of  applicability  of  our  sensitivity 
analysis.  Since  some  of  our  input  parameters  have  wide  limits  of  uncertainty, 
the  sensitivity  of  the  system  may  change  considerably  as  our  parameter  value 
is  changed.  It  is  therefore  useful  to  have  some  idea  of  the  region  in  which 
the  sensitivity  analysis  is  valid. 

In  this  paper  the  results  are  given  for  a  sensitivity  analysis  performed 
on  a  set  of  H2  -  O2  -  N2  flames.  The  input  parameters  considered  variable  are 
the  transport  coefficients  and  the  rate  coefficients.  Since  the  thermodynamic 
properties  for  this  flame  are  well  known,  our  sensitivity  analysis  has  not 
been  extended  to  cover  these  input  parameters.  Similarly,  because  the  ratio 
between  the  forward  and  the  backward  rate  coefficients  of  a  reaction  is 
determined  by  the  thermodynamic  properties,  the  effect  of  varying  this  ratio 
is  not  considered  either.  The  appendix  outlines  the  sources  of  input 
parameters  used. 

Section  II  gives  the  governing  equations,  both  for  the  basic  model  and 
for  the  sensitivity  coefficients.  Section  III  reports  the  results  for  the  two 
special  cases  that  can  be  solved  analytically.  In  Section  IV,  we  discuss  a 
numerical  procedure.  Section  V  gives  the  numerical  results  for  our  test  set 
of  flames.  In  Sections  VI  and  VII  we  develop  the  second  order  sensitivity 
analysis,  and  in  Section  VIII  we  discuss  the  accuracy  of  extrapolating  the 
results  when  the  parameters  vary  over  a  large  range. 

II.  GOVERNING  EQUATIONS 

We  are  interested  in  the  equations  that  describe  a  one-dimensional, 
laminar,  preraixed  flame  propagating  in  an  unbounded  ideal  gas.  Since 
performing  the  sensitivity  analysis  numerically  is  expensive,  it  is  useful  to 
simplify  the  equations  as  much  as  possible.  A  high  degree  of  accuracy  is  not 
required,  since  we  are  interested  in  trends. 


o 


In  a  previous  paper^,  we  considered  various  levels  of  approximation  to 
the  multicomponent  transport  coefficients.  A  relatively  simple  approximation 
(Method  V),  gave  very  accurate  results  and  will  be  used  here.  The 
corresponding  equations  are  the  following. 


Continuity  for  species  1  through  N~1  is  expressed  in  the  form 

aY.  aY.  .  32Y,  R.M, 

+  m_  ttt-  =  P  D^_  — —  +  — —  i  =  1,  N-l,  (1) 


3t 


o  34- 


im  34, 2 


and  continuity  of  the  Nth  species  in  the  form 


N-l 

Y  ■  1  -  E  Y 

N  i-1  i 


(2) 


Here  t  is  the  time,  Y.  the  mass  fraction  of  the  ith  species,  raQ  the  mass  flux 
through  the  cold  boundary;  ^  a  transformed  space  coordinate,  P  the  density  of 
the  mixture,  Djm  the  diffusion  coefficient  of  the  ith  species  into  the 
mixture,  the  molecular  weight  of  species  i  and  the  net  rate  of 
production  of  species  i  due  to  chemical  reactions.  The  coordinate  is 
related  to  the  spatial  coodinate  x  by 

X  /  / 

4-  =  /  P  (x  )  dx  .  (3) 

o 

Conservation  of  energy  is  expressed  by 


3T  3T  pX  32T 

3t  o  34-  c  a.2 

P  3+ 


R.M.h.  , 
i  1  i 


(4) 


where  T  is  the  temperature  of  the  mixture,  X  the  thermal  conductivity  of  the 
mixture,  c  the  specific  heat  of  the  mixture  and  h^  the  specific  enthalpy  of 
species  i,  given  by 


h. 

1 


h.°  +  c  (T  -  T  ) 
1  p  o 


(5) 


where  Tq  is  a  fixed  reference  temperature  and  h^°  is  the  enthalpy  at  that 
temperature.  The  boundary  conditions  are 


T  =  T  and  Y,  =  Y.  ,  i  =  1,2,  ...  N-l 
u  i  iu 


(6) 


2r.P.  Coffee  and  J.M.  Heimerl ,  "Transport  Algorithms  for  Premixed ,  Laminar 
Steady  State  Flames ,  "  Combustion  and  Ftamej  Vol.  43,  pp.  273-289,  1981. 
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(7) 


at  x  =  ^  = 


9T 

3^ 


-  00  and 


3  Y. 


i 


1,2, . ..N-l 


at  x  8  \[i=4*00.  The  above  system  of  equations  is  closed  by  relations  that 

express  the  production  rates  in  terms  of  the  temperature,  mass  fractions, 
and  a  set  of  chemical  parameters. 

In  the  above  equations,  the  effects  of  radiation,  viscosity,  and  body 
forces  are  ignored.  Since  the  burning  velocity  is  small  compared  with  the 
local  speed  of  sound,  the  pressure  is  taken  to  be  constant.  Besides  these 
standard  assumptions,  we  have  also  assumed  that  the  thermal  diffusion  is 
negligibly  small,  that  Fick/s  law  holds  for  the  diffusion  velocities,  and  that 
the  quantities  p^D^m,  pX  ,  and  c^  are  constant.  These  constants  are  chosen  a 
priori . 

Our  interest  is  in  the  steady  state  solution.  The  numerical  procedure  is 
to  start  with  arbitrary  profiles  Y^(^)  and  T(^)  and  integrate  in  time  until 
the  steady  state  solution  is  achieved.  As  the  integration  proceeds,  the  mass 
flux  m  through  the  origin  is  iteratively  adjusted  to  equal  the  mass  flux 
through  the  flame,  so  that  at  steady  state  the  flame  is  motionless  with 
respect  to  the  coordinate  system  used.  The  burning  velocity,  S,  of  the  flame 
is  defined  as  the  velocity  of  the  flame  relative  to  the  fluid  at  rest,  i.e., 
at  x  =  -  00 .  Since  in  our  coordinate  system  the  flame  is  at  rest,  S  =  v(-  00 )  = 
mn/  Pu,  where  Pu  is  the  density  of  the  unburned  mixture. 

The  solution  of  the  Eqs .  (1)  through  (7)  establishes  benchmark  values  of 
the  burning  velocity  and  the  temperature  and  species  profiles  for  a  given  set 
of  input  parameters.  Details  of  the  numerical  procedure  are  in  References  2 
and  3 . 


The  production  rate  functions,  R^(T),  are  defined  by 

NR  NR  N 

E 


NR  NR  N  ,  N  „ 

‘  r“l  hr  -  r=l  <h,r  '  V  i.r'>  [kfr  >1  <°J> * *  3’r  "Sr  jSl^5  j>rl 


(8) 


where  v" •  and  v ' ^  are  the  number  of  molecules  of  species  i  entering  into 
reaction  r  as  product  and  reactant,  respectively,  Cj  is  the  concentration  of 
the  jth  species  and  k^  and  k^  are  the  forward  and  reverse  rate  coefficients, 
respectively .  The  rate  coefficients  depend  on  the  temperature  through  the 


%\P.  Coffee  and  J.M.  Heimevl,  "A  Method  for  Computing  the  Flame  Speed  for  a 
Laminar ,  Premixed,  One  Dimensional  Flame, "  BRL  Technical  Report,  ARBRL-TR- 

02212,  Jan.  1980,  (AD#A082803) . 

3T.P.  Coffee,  "A  Computer  Code  for  the  Solution  of  the  Equations  Governing  a 
Laminar,  Premixed,  One  Dimensional  Flame,  BRL  Memorandum  Report,  ARBRL-MR- 
03165 ,  April  1982,  ( AD#A114041 )  . 
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Arrhenius  relations 


kfr  =  Afr  T  Bfr  exp(-Efr/RT) 


and 


kbr  =  Abr  T  Bbr  exP(”Ebr/RT) ,  (9) 

where  R  is  the  gas  constant,  E  is  the  activation  energy  and  A  is  the  frequency 
factor  for  the  reaction.  In  Eg.  (9)  the  parameters  Bfr,  Bbr,  Efr  and  E^r  will 
be  held  constant  and  sensitivity  analysis  will  be  restricted  to  variations  in 
Afr«  Since  thermodynamic  quantities  are  held  constant  the  ratio  is 

also  constant.  We  will  also  investigate  the  effects  of  variations  of  the 
parameters  p2D^m,  i  =  1,  . ..,  N  -  1  and  pA  in  Eqs .  (1)  and  (4).  The 
remaining  parameters  Cp,  and  h^°  are  held  fixed.  Thus  we  have  a  total  of 
NT  =  N  +  NR  variable  parameters  .  See  Table  1  . 


*  TABLE  1  .  PARAMETERS  DEFINING  THE  FLAME. 


Par a me ter 
Mixture 

Species  (i  =  1,  • • • ,  N) 

Reaction 
( r=l , . • • , NR) 


Constant 

cp 

h°  , 

1  1 

Efr'  Bfr 
^r'^r 


Variable 


pA  x  a 


NT 


p  D.  x  ot 

lm  NR  +  i 

x  ar 


Abr  x  ar 


Hie  variation  is  expressed  by  attaching  factors  a.,  i 
the  benchmark  values  of  the  variable  parameters.  Then  feqs. 


=  1,  . . NT  to 

(1)  and  (4)  become 


.Y.  Y. 

3  l  3 _ l 

3t  +  mo  3\|> 


2 

=  a  .  p  D. 
NR-fi  lm 


32y.  nr 

- 1+  E  ,  a  R.  M. /p ,  i=l, 2. 

. , 2  r=l  r  lr  1 

3^ 


.  ,  N— 1 


(10) 


and 


3  T  3  T  aNT  PX  3  2T 

3 1  +  mo  3ty  c  »  .  2 

p  3^ 


N  NR 

—  Z  ,  E  a  R.  M.  h.  .  (11) 

c  1=1  r=l  r  lr  l  l 
P 


When  a  =  . .  .  =  a  NT=1/  Bqs .  (10)  and  (11)  are  just  the  benchmark  Eqs.  (1) 
and  (4)1  Changes  in  a  x  (i=l,  ...  NR)  correspond  to  changes  in  the  rate 
coefficients  kfi  and  k^  ?  changes  in  a  (i  =  NR+1,  ...  NT)  correspond  to 

2 

changes  in  the  transport  coefficients  p  D  .  and  pA  .  The  basic  problem  is 
to  determine  how  changes  in  the  affect  tRe  output  functions  S,  and  T. 
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III.  ANALYTIC  RESULTS 


Two  special  cases  can  be  solved  analytically.  Consider  the  case 

a,  =  ...  aTO  =  1,  a  =  ...  =  a  =  a  Then  Eqs.  (10)  and  (11)  can  be 

1 .  NR  NR+1  NT 

written 


3Yi 

7T  +  “o  ~4 


.  32Y.  R_,M, 

l  ,  i  i 

ap  D.  - »  +  — — 

im  2  p 


(12) 


and 


3T  ,  3  T 

TT  +  mo  W 


2 

otpA  8  T 

c  a  .  2 

P  3 'I' 


1 

Pc 

P 


U 


R.  M. 
l  i 


(13) 


Apply  the  coordinate  transformation  =  ty/fuT . 

become : 


~Y. 

8  l 
8 1 


m 


+ 


/a 


iii 

r  9*' 


2  ,  Vl 

#  ^  +  T- 


Then  the  above  equations 


(14) 


and 


3T  m°  3T_  =  pX_  3^T 

8 1  7a  dty'  c  ^ .  ,  2 

P  H 


Pc 


N 

iz=i 


\  Mi  V 


(15) 


Eqs,  (14)  and  (15)  are  formally  identical  to  the  benchmark  Eqs.  (1)  and  (4) 
except  that  the  mass  flux  factor  mQ  in  these  equations  is  replaced  by 
m  /  /a.  Thus,  multiplying  all  the  transport  coefficients  by  a  uniform  factor 
a  causes  the  mass  flux  (and  hence  the  burning  velocity)  to  change  by  a  factor 
/a.  Similarly,  the  profiles  in  space  are  identical  to  the  benchmark 
profiles  in  ^  space.  That  is,  the  profile^  are  expanded  (a  >  1)  or 
contracted  (a  <  1)  uniformly  by  a  factor  /a  in  the  original  space.  This 
fact  has  been  noted  by  Dixon-Lewis^ . 


A  similar  result  follows  if  all  the  rate  coeffients  are  multiplied  by  a 
uniform  factor  a,  that  is,  =  • • •  =  a >aNR+l  =  •••  =  aNT  ~  Then 

using  the  coordinate  transformation  -  ^/cx  and  dividing  by  a,  one  obtains 

(16) 


n  ~Y,  m  0Y-r 

1  3  i  +  _o  8  i 


a  8t 


/a  8ij> 


2 

2.  9  h  , 

=  P  Dim  7^2  + 


R .  M , 
i  i 


8i() 


4G.  Dixon-Lewis,  "Flame  Structure  and  Flame  Reaction  Kinetics.  I.  Solution  of 
Conservation  Equations  and  Application  to  Rich  Hydrogen-Oxygen  Flames ,  Proc . 
Rojy .  Soc_.  ,  London  A,  Vol .  298 ,  pp .  495-513,  1967. 
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and 


l  3T  ,  mo  3T 
a  3t  /—  S'P ' 


pX_  j^r 

c  a  ,  rl 
P  9^ 


1 

PCP 


u  RiMihi 


(17) 


For  the  steady  state  solution,  multiplying  all  the  raj^es  by  a  uniform  factor 
a  causes  the  burning  velocity  to  change  by  a  factor  /a  while  the  profiles  are 
contracted  (a  >  1)  or  expanded  (a  <  1)  in  the  same  proportion. 


So  the  overall  effect  of  increasing  the  transport  coefficients  is  to 
spread  the  flame  front,  as  energy  and  species  diffuse  more  rapidly  away  from 
the  flame  front.  The  effect  of  increasing  all  the  rates  is  to  contract  the 
flame  front,  as  the  overall  combustion  occurs  more  rapidly. 


These  effects  are  independent.  That  is,  if  all  the  transport 
coefficients  and  all  the  rate  coefficients  are  multiplied  by  the  same  factor 
a  ,  then  the  flame  speed  S  will  change  by  a  factor  a,  while  the  species  and 
temperature  profiles  will  be  identical  with  the  solutions  of  the  benchmark 
equations . 


IV.  NUMERICAL  PROCEDURES 

Several  approaches  have  been  used  to  perform  sensitivity  analyses,  mostly 
for  systems  of  ordinary  differential  equations  in  time.  These  are  either 
local  (valid  for  input  parameters  near  the  benchmark  values)  or  global  (valid 
over  a  specified  range).  Since  some  rate  coefficients  have  large  regions  of 
uncertainty,  a  global  procedure  would  be  preferred.  The  simplest  procedure  is 
to  solve  the  equations  for  different  values  of  the  input  parameters  over  their 
range  of  uncertainty.  For  large  numbers  of  parameters  whose  regions  of 

t .  A  more 
,  the  Fourier 

Amplitude  Sensitivity  Test  (FAST).  In  this  procedure  the  input  parameters  are 
varied  simultaneously,  and  the  results  are  Fourier  analyzed.  However,  as  the 
number  of  input  parameters  increases,  the  number  of  computations  required 
increases  rapidly  and  the  cost  for  our  system  of  PDE's  is  still  prohibitive. 

Therefore  we  use  a  local  method  based  on  a  Taylor  series  expansion  with 
respect  to  variations  of  the  parameters  a  . .  We  introduce  the  notation 


uncertainty  are  also  large,  this  simple  procedure  is  inefficien 
efficient  procedure  has  been  devised  by  Shuler  and  coworkers^ 


^R.I.  Cukier,  C.M.,  Fortuin,  K.E .  Shuler,  A.G .  Petschek,  and  J.H .  Schaibley, 

" Study  of  the  Sensitivity  of  Coupled  Reaction  Systems  to  Uncertainties  in  Rate 
Coefficients  I.  Theory ,'  &U  Chem*  .  Vol .  59,  pp.  3873-3878,  1973. 

°J.H.  Schaibly  and  K.E.  Shuler,  "Study  of  the  Sensitivity  of  Coupled  Reaction 
Systems  to  Uncertainties  in  Rate  Coefficients  II.  Applications,  d.  Chenu 
Ehy£.  •  •  Vol.  59,  pp.  3879-3888,  1973. 

?R.I.  Cukier,  J.H.  Sahaibly,  and  K.E.  Shuler,  "Study  of  the  Sensitivity  of 
Coupled  Reaction  Systems  to  Uncertainties  in  Rate  Coefficients  III.  Analysis 
of  the  Approximations >  tL  Chem.  Vol.  63,  pp.  1140-1149,  1975. 


li 


F^  =  9F/9a  .  evaluated  =  ...  =  a  =1  (F 

flame  speed,  S,  to  first  order  and  find 
NT 


SN“ 


SE+  5-1  S 


-  S,  or  T) .  We  expand  the 


(18) 


where  is  the  benchmark  flame  speed  and  SN  is  the  flame  speed  for  other 
choices  of  the  a^.  Similar  expansions  are  used  for  and  T. 

The  partial  derivatives  ,  Y^  and  T^ ,  i.e.  the  "sensitivity 
coefficients",  are  computed  numerically.  To  that  end  we  take  the  partial 
derivatives  of  Eqs.  (10)  and  (11)  with  respect  to  ot  .  and  let  the  ot's  approach 
one.  The  resulting  equations  are  ^ 


3Y.J  .  3Y3  3Y.j  „  32Y.  82Y  j 

~~  +  m  3  +  m  =  6  .  .  P  D,  - K  +  P  D  - ,  + 

at  o  9\J)  o  3\J)  NR4-1 ,  j  im  ^.,2  im  ^.,2 


9iJ) ‘ 


R,  M,  M.  R, 

F  [6  .  _i!U_  +  (J-Ji  )  3] 

r=l  r,j  p  p 


9\p 
(19) 


and 


3T^ 


j  3T  .  9  T 


^ —  4-  m  ttt-  4-  m 

3 1  o  3i|i 


-  pX  3^T  pX_  32Tj 

o  3*  =  NTJ  cp  ^2  cp  ^  2 


.  N  NR  6  .  R.  M.  h. 

1_  y  v  r  r>J  lr  1  1 

c  1=1  r=l  1  p 

P 


Rir  Mi  hi  3  Rir  Mi  ^ 

+  -p  —  +  V 


(20) 


where  from  Eq .  (5)  we  see  that  h^  =  c  T 3 .  Introduction  of  the  notation 


u±  =  Yi,  i  =  1,  ...N, 
and  uN+1  =  T, 


(21) 


allows  us  to  write  the  identity 

R.  M.  .  N+l  R  M 

(  ir  X~)  3  =  E .  [3  (  lr  ■  ■  — ) / 3 u  ] u  3  .  (22) 

v  p  '  m=l  p  m  m 

Equations  (19)  through  (22)  constitute  a  system  of  differential  equations 
for  the  functions  u  This  system  is  similar  to  that  given  by  Eqs.  (1) 

through  (7)  and  is  solved  using  the  same  numerical  relaxation  technique 
outlined  in  Section  II. 

This  numerical  solution  is  called  the  Direct  Method  and  is  accomplished 
as  follows.  First  the  benchmark  Eqs.  (1)  and  (4)  are  solved.  The  benchmark 
values  for  mQ ,  the  values  for  the  species  and  temperature,  and  their  first  and 
second  derivatives  are  stored  in  a  file.  (We  use  a  finite  element  code  in 
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which  the  solutions  are  represented  by  piecewise  polynomials.  Therefore  the 
derivatives  are  available.  Also  the  code  uses  a  collocation  method  so  that 
only  the  values  at  the  collocation  points  need  to  be  saved.)  The  terms 
(RirMf/p)  and  9  (R^frh/p  )/9 u^  are  computed  in  subroutines  and  the  benchmark 
values  of  these  expressions  are  also  saved.  Finally,  Eqs.  (19)  and  (20)  are 
solved  for  each  selected  set  of  a.,  providing  us  with  the  corresponding  set  of 
mQ^  and  u^  .  ^ 

Another  method  for  finding  the  linear  sensitivity  coefficients  is  the 
Green’s  function  method.^  For  ordinary  differential  equations  this  involves 
solving  a  set  of  N  equations  to  determine  the  Green’s  functions,  instead  of  a 
set  of  NT  equations.  Then  the  sensitivity  coefficients  are  found  by  a  set  of 
quadratures.  Higher  order  sensitivity  coefficients  can  also  be  determined  by 
the  same  method. 

Q 

In  a  recent  paper  ,  the  theory  was  extended  to  partial  differential 
equations.  In  that  case  the  Green’s  functions  are  represented  as  a  linear 
combination  of  independent  solutions  of  the  homogeneous  equations  based  on  our 
Eqs.  (19)  and  (20).  Each  independent  solution  is  found  by  solving  the 
homogeneous  equations  with  a  different  set  of  boundary  conditions.  We 
implemented  this  procedure,  but  the  solutions  obtained  were  not  numerically 
independent.  This  is  due  to  the  fact  that  the  steady  state  solutions  of  the 
flame  equations  are  almost  independent  of  the  downstream  (or  hot)  boundary 
conditions.  So  this  procedure  cannot  be  used  for  the  steady  state  flame 
equations . 

Once  the  sensitivity  coefficients  have  been  obtained,  the  usual  way  to 
determine  the  values  of  the  output  functions  that  correspond  to  changes  in  the 
input  parameters  is  by  the  linear  terms  of  a  Taylor  series,  such  as  Eq. 

(18).  However,  the  results  of  Section  III  indicate  that  at  least  for  two 
special  cases,  the  extrapolation  S^  =  Sg  a’5  is  exact.  This  suggests  that  a 
better  approximation  than  Eq.  (18)  may  be  the  formula 

NT  s  j 

SX7  —  S  n  ,  o.  E  .  (23) 

N  B  j=l  j 

where  the  S£J  are  appropriate  sensitivity  exponents.  Taking  logarithms,  Eq. 
(23)  becomes 


8J.T.  Huang,  E.P.  Dougherty,  S .  Rabitz  and  H .  Babitz,  " The  Green1 8  Function 
Method  of  Sensitivity  Analysis  in  Chemical  Kinetics J.  Chem .  Phys.,  Vol. 
69,  pp.  5180-5191,  1978. 


9M.  Demiralp  and  H .  Babitz,  " Chemical  Kinetic  Functional  Sensitivity 
Analysis:  Elementary  Sensitivities*  J .  Chem.  Phys.%  Vol.  74,  pp.  3362-3375 . 
1981. 
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In 


SN 


“  In 


NT 

+ 1-1 


lna  . 
J 


(24) 


Eq.  (24)  represents  the  linear  Taylor  expansion  of  In  (SN)  in  terms  of  In  a^ . 
This  form  permits  some  checks  on  the  numerical  solution.  Eq.  (23)  becomes 


exact  as  the  a/s  approach  one.  But  if  we  consider  th^  case  where 
“l-  •••  -“m°r  “NR+l  ■  "•  -°OT>  then  SN-  SB  •  E’uatln 


we  If ind 


Equating  exponents, 


NR 

i-i 


J 


E 


NT 

j=NR+l 


E 


J  =  0.5, 


(25) 


The  corresponding  relation  for  the  profile  sensitivity  coefficients  is 
NT 


5-i 


UiEJ  "° 


(26) 


for  any  value  of 


The  relation  between  the  sensitivity  coefficients  and  is 

straightforward.  Taking  the  derivative  of  Eq.  (24)  with  respect  to  am 
results  in 


8a 

m 


and  in  the  limit  a  +  1 , 

ra 


(27) 


S  /  Sr>  . 


(28) 


The  corresponding  formula  for  the  u 


j 

iE 


is 


u 


j 

iE 


(29) 


The  logarithmic  sensitivity  form  (24)  has  been  used  before  .  In 
particular,  it  has  been  recommended  for  the  modelling  of  chemical  reactions 
without  transport,  where  logarithmic  measures  of  species  concentrations  as  a 


•in  * 

10P.M.  Frank,  "Introduction  to  System  Sensitivity  Theory,  Academic  Press, 
1978. 
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function  of  time  are  appropriate11.  We  will  show  that  it  is  also  useful  for 
steady  state  flames. 

V.  NUMERICAL  RESULTS  AND  DISCUSSION 

The  above  procedure  was  carried  out  for  a  set  of  four  ~  ~  N2 

flames.  The  chemistry  scheme  used  is  due  to  Dixon-Lewis1  and  consists  of 
eight  species  and  seventeen  reactions.  The  input  parameters  are  given  in 
Appendix  A. 

The  initial  mole  fractions  and  temperature  for  the  flames  are  given  in 
Table  2.  The  benchmark  flame  speeds  are  also  given.  These  values  are 
slightly  different  from  the  flame  speeds  calculated  using  a  more  accurate 
transport  algorithm.  The  total  pressure  is  fixed  at  one  atmosphere  for  all 
the  flames. 

Flame  A  is  a  very  fast,  stoichiometric  flame.  Flame  B  is  a  hydrogen  rich 
flame  with  a  more  moderate  flame  speed.  Flame  C  is  an  oxygen  rich  flame.  The 
initial  mole  fractions  were  chosen  so  that  the  flame  speed  is  similar  to  that 
of  flame  B.  Flame  D  is  a  very  slow,  hydrogen  rich  flame.  This  flame  has  been 
studied  extensively  by  Dixon- Lewis1  . 

TABLE  2.  MOLE  FRACTIONS,  INITIAL  TEMPERATURE  AND  BENCHMARK  FLAME 
SPEEDS  FOR  THE  FOUR  FLAMES  STUDIED. 


FLAME 

XH2 

CM 

o 

X  ! 

1 

CM 

VK) 

Sg(cm— 

A 

.6667 

.3333 

.0000 

298 

858.0 

B 

.5000 

.1050 

.3950 

298 

251.9 

C 

.3000 

.5600 

.1400 

298 

285.7 

D 

.1883 

.0460 

.7657 

336 

11.8 

21W.C.  Gardiner,  Jr.,  "The  pC,  pR,  pP,  pM,  and  pS  Method  of  Formulating  the 
Results  of  Computer  Modeling  Studies  of  Chemical  Reactions,  J.  Phys .  Cherry, , 
Vol  .  81,  pp.  2367-2871,  1977. 

Dixon-Lewis,  " Kinetic  Mechanism,  Structure  and  Properties  of  Premixed 
Flames  in  Hydrogen-Oxygen-Nitrogen  Mixtures,"  Proa .  R.  Soc .  London  A,  Vol. 
292,  pp.  45-99,  1979. 

Dixon-Lewis,  M.M.  Sutton,  and  A.  Williams,  "Flame  Structure  and  Flame 
Reaction  Kinetics  IV.  Experimental  Investigations  of  a  Fuel-Rich  Hydrogen  + 
Oxygen  +  Nitrogen  Flame  at  Atmospheric  Pressure, 

Proc.  Ron.  Soc.  London  A.  .  Vol.  312,  pp.  227-234,  1970. 
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Table  3  gives  the  sensitivity  exponents  of  the  transport  parameters 
and  Table  4  gives  them  for  the  rate  coefficients.  Where  possible  we  have 
taken  the  Nth  specie^  of  Eq .  (2)  to  be  the  diluent,  Thus,  for  flames  B,  C 

and  D  an  entry  for  p  does  not  appear.  (See  Table  1).  Since  flame  A  has 

no  diluent  and  is  described  by  one  less  species,  we  have  taken  the  Nth  species 
to  be  H2O.  Consequently,  no  entry  appears  in  Table  3  for  ^0,  flame  A. 
Likewise  in  Table  4  reaction  12  is  zero  and  no  entry  appears. 

Comparing  the  E  S  ^  from  the  numerical  values  of  Tables  3  and  4  with  the 
analytic  result  given  by  Eq .  (25),  we  find  excellant  agreement. 

The  tables  also  give  the  sum  of  the  absolute  values  of  the  S^ •  We  will 
refer  to  this  number  as  the  overall  sensitivity  (for  transport  or 
chemistry).  It  will  be  shown  that  this  number  provides  a  rough  estimate  of 
the  complexity  of  the  flame  and  its  sensitivity  to  changes  in  the  input 
parameters • 

Consider  first  the  transport  parameters.  A  flame  front  can  propagate  by 
two  processes:  diffusion  of  radicals  (H,0H,0,H02)  ahead  of  the  flame  or 
conduction  of  heat  ahead  of  the  flame.  Among  the  radicals,  the  flames  are 
most  sensitive  to  the  diffusion  of  H.  This  is  reasonable  since  H  is  light  and 
diffuses  easily.  The  faster  flames  are  more  sensitive  to  the  diffusion  of 

TABLE  3.  FLAME  SPEED  SENSITIVITY  EXPONENTS,  S£j  , 

FOR  THE  TRANSPORT  PARAMETERS. 


Parameter 

,A 

B 

C_ 

D 

H 

.3545 

.3785 

.2018 

.0123 

p2  D2m 

OH 

.0250 

.0012 

.0647 

-.0013 

p2°3m 

0 

.0339 

-.0011 

.0936 

-.0033 

p\„ 

ho2 

.0001 

-.0000 

-.0008 

.0001 

PS„ 

h2 

-.0069 

.0079 

-.2689 

-.0025 

p2°6ra 

°2 

-.0237 

-.1207 

-.0159 

-.4176 

2„ 

P  D7m 

h2o 

* 

.0162 

.0232 

-.0406 

pX 

.1179 

.2191 

.4039 

.9521 

NT 

e  sEi 

j=NR+l 

.5009 

.5010 

.5017 

.4992 

NT 

£  lSEJ| 

j=NR+l 

.5620 

.7446 

1.0727 

1.4298 

*no  entry,  see  text. 
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TABLE  4.  FLAME  SPEED  SENSITIVITY  EXPONENTS,  , 

FOR  THE  RATE  COEFFICIENTS 

Race  Coefficient  Reaction  ABC  D 


1 

0H+H2 

-  h2o+h 

.1346 

.0744 

.3140 

.0576 

2 

H+02  **  OH+O 

.1041 

.2911 

.0400 

1.2427 

3 

0+H2  ~  OH+H 

.0907 

.0468 

.1718 

.0347 

4 

H+02+M' 

-  ho2+m' 

.1879 

.1838 

.0301 

-.5312 

3 

h+ho2 

-  OH+OH 

.0503 

.1099 

.0997 

.7492 

6 

h+ho2 

• 

-  o+h2o 

-.0053 

-.0054 

-.0111 

.0389 

7 

h+ho2 

-  h2+o2 

-.0399 

-.0952 

-.0666 

-.7555 

8 

oh+ho2 

-  h2o+o2 

-.0056 

-.0092 

-.0208 

-.0187 

9 

OHIO  2 

-  oh+o2 

-.0005 

-.0011 

-.0006 

-.0125 

10 

ohio2 

-  oh+o2 

-.0001 

-.0002 

-.0000 

-.0022 

11 

h+h+h2 

-  h2+h2 

-.0003 

-.0142 

-.0002 

-.0222 

12 

h+ii+n2 

-  h2+n2 

* 

-.0111 

-.0001 

-.1146 

13 

h+h+o2 

"  h2+02 

-.0001 

-.0009 

-.0005 

-.0007 

14 

h+h+h2o 

-  h2+h2o 

-.0014 

-.0419 

-.0014 

-.1339 

15 

H+OH+M" 

-  h2o+-m" 

-.0047 

-.0213 

-.0112 

-.0296 

16 

H+O+M" 

-  OH+M" 

.0001 

-.0007 

-.0008 

-.0015 

17 

OH+OH 

-  o+-h2o 

-.0106 

-.0060 

-.0440 

.0001 

NR 

z  sEJ 
j=l 

.4991 

.4990 

.4983 

.5008 

NR 

z  lsPj  ! 

.6362 

.9131 

.8128 

3.7458 

j=l 

*  no  entry,  see  text. 
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radicals,  since  their  flame  fronts  have  larger  radical  concentrations  and 
steeper  gradients  than  the  slower  flames.  Increasing  the  thermal  conductivity 
increases  the  flame  speed.  This  effect  progressively  becomes  more  important 
the  slower  the  flame  because  diffusion  becomes  less  important. 

The  diffusion  of  R2  and  O2  tends  to  slow  down  the  flame,  since  they  tend 
to  diffuse  away  from  the  flame  front.  The  hydrogen  rich  flames,  B  and  D,  are 
sensitive  to  the  diffusion  of  oxygen  since  oxygen  is  in  short  supply,  while 
the  oxygen  rich  flame,  C,  is  sensitive  to  the  diffusion  of  hydrogen. 

The  results  for  the  chemistry  are  consistent  with  the  prevailing  view  of 
the  mechanism.  Most  of  the  product  1^0  is  created  by  reaction  1 
[OH+H2  ~  R2O+R]  and  increases  in  this  rate  will  increase  the  flame  speed. 

Note  that  flame  C  is  especially  sensitive  to  this  rate.  Since  flame  C  is 
deficient  in  R2,  reaction  1  is  slower  relative  to  the  other  flames.  A  slower 
reaction  can  be  a  bottleneck  for  an  entire  process  with  the  consequence  that  a 
flame  can  become  quite  sensitive  to  changes  in  the  bottleneck  reaction  rate. 

The  OR  required  in  reaction  1  is  primarily  produced  by  the  chain¬ 
branching  reactions  2  [R+O2  ~  0R+0]  and  3  [OHI2  **  OR+H] •  The  hydrogen  rich 
flames  B  and  D  are  more  sensitive  to  reaction  2  (shortage  of  O2)  and  the 
oxygen  rich  flame  C  is  more  sensitive  to  reaction  3  (shortage  of  R2) •  F°r  the 
stoichiometric  flame  A,  the  sensitivities  are  comparable. 

OR  can  also  be  created  by  the  reaction  pathway  4[R+02+M  ~  RC^+M]  and 
5[R+R02  ~  0R+0R] .  Reaction  5  is  in  competition  with  reaction 
7 [ R+R02  ~  For  flames  A,B ,  and  C,  increasing  rate  coefficients  4  and  5 

will  increase  the  flame  speed,  while  increasing  rate  coefficient  7  will 
decrease  the  flame  speed.  The  other  reactions  are  relatively  unimportant  for 
these  flames. 


Because  of  the  fact  that  the  concentration  of  R  is  very  small,  flame  D 
exhibits  a  different  behavior.  If  rate  coefficient  4  is  increased,  the  flame 
speed  decreases.  That  is,  the  shortage  of  R  is  aggravated,  which  slows  down 
the  total  network.  For  the  same  reason,  flame  D  is  very  sensitive  to  the 
recombination  reactions  (11,12,  and  14). 

Flame  D  is  in  fact  very  close  to  the  point  of  extinction.  A  relatively 
small  change  in  an  input  parameter  (increasing  rate  coefficients  4  or  7 ; 
decreasing  rate  coefficients  2  or  5)  can  extinguish  the  flame.  Conversely,  a 
relatively  small  change  in  the  opposite  direction  can  lead  to  a  large  increase 
in  the  flame  speed.  This  complex  balance  between  several  competing  reactions 

NR 

is  reflected  in  the  fact  that  the  overall  sensitivity  £  | SEJ |  is  large. 

j  =  l 


We  are  also  interested  in  the  sensitivity  of  the  species  and  temperature 
profiles.  In  this  case  there  are  some  problems  in  presenting  the  results  in 
terms  of  the  sensitivity  exponents  given  by  Eq.  (29).  Near  the  upstream  (or 
cold)  boundary,  some  of  the  mass  fractions  u^  approach  zero.  Thus  in  this 


region  small  numerical  errors  in  either. u^ 
As  a  heuristic  correction  we  set  u“  1= 


iE 


iE 


or  u.  lead  to  large  errors  in 
0  ifu  <  ui^(max)  x  10”^  where 
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Ui^(max)  is  the  maximum  value  of  the  benchmark  profile.  Another  difficulty 
occurs  just  ahead  of  the  flame  front.  Our  analytic  results  indicate  that  a 
flame  front  may  expand  due  to  changes  in  the  input  parameters.  In  that  case, 
a  small  value  of  u^  just  ahead  of  the  flame  front  can  correspond  to, a 
relatively  large  sensitivity  coefficient  u^ ,  and  therefore  the  u.^  *-s 

large.  This  large  value  can  determine  the  scale  of  a  plot  of  the  logarithmic 
sensitivity  profiles,  with  the  result  that  the  details  of  these  profiles  in 
the  region  of  the  flame  front  may  be  lost.  Therefore,  it  is  more  practical  to 
graph  the  sensitivity  coefficients  u^  rather  than  the  u 

For  ease  of  comparison,  we  normalize  the  profiles  by  dividing  the  u^  by  the 
u  (max) •  (All  calculations  are  carried  out  in  ^  -  space  and  ,  for  easier 
interpretation,  are  converted  to  x-space  before  graphing.) 


Figures  1  and  2  show  the  benchmark  mass  fraction  profiles  for  OH 
(radical)  and  H2  (major  species),  respectively.  We  have  taken  these  two 
species  and  flame  B  as  representative  examples.  The  normalized  sensitivity 
coefficients  profiles  for  p  D  ,  p  T>  ,  p^D  ~  and  pX  are  given  in  Figures 
and  4.  The  profiles  for  reaction  rate  coerr icients  1,2, 4, 5,  and  7  are  shown 
in  Figures  5  and  6.  All  other  transport  and  rate  coefficient  sensitivity 
profiles  for  flame  B  are  small.  The  profiles  for  OH  are  larger,  because  of 
the  normalization.  A  small  relative  change  in  a  major  species  concentration 
can  lead  to  large  relative  changes  in  the  radical  concentrations,  because 
the  radical  concentrations  are  much  smaller. 
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In  general  the  diffusion  of  H  is  the  most  important  transport 
parameter.  The  sensitivity  coefficients  in  Figure  3  show  that  increasing 
p  leads  to  a  sharp  increase  of  the  first  OH  bump  in  the  early  part  of  the 

flame  of  Figure  1.  We  see  in  Figure  4  that  the  diffusion  of  H2  is  important 
in  determining  the  H2  profile,  but  it  has  a  minor  effect  on  all  the  other 
profiles . 

The  chemistry  profiles  are  more  complicated,  since  different  reactions 
can  be  important  at  different  locations  in  the  flame  front.  In  Figures  5  and 
6  we  see  that  rate  coefficient  2  (F+H2  ~0H+0)  is  the  most  important.  The 
complicated  appearance  of  the  sensitivity  coefficients  for  the  OH  profiles  in 
Figure  5  is  due  to  competition  between  the  chain  branching  rates  and  the  HO2 
rates  in  the  early  part  of  the  flame  front. 

The  results  for  flame  C  (not  shown)  are  similar  to  those  given  above, 
except  that  rate  coefficient  3(0+^  ~0H+H)  becomes  important. 

For  the  fast  flame  A,  the  radical  concentrations  are  almost  as  large  as 
the  major  species  concentrations  (YQ^^max^=  .1104;  YH2^maX  -  .1111).  As  a 
result,  the  normalized  OH  sensitivity  coefficients  are  about  the  same  size  as 
the  H2  coefficients.  The  profiles  are  relatively  insensitive  to  changes  in 
the  input  parameters.  In  general,  the  same  set  of  input  parameters  as  before 
are  important  for  flame  A.  The  exception  is  in  the  post  flame  region  where, 
because  of  the  large  radical  concentrations,  rate  coefficients  13  and  14 
become  important • 

For  the  slow  flame  D,  the  radical  concentrations  are  very  small  (YQ^max) 
=  1.09  x  10  0,  so  that  the  normalized  sensitivity  coefficients  are  quite 
large  (on  the  order  of  1.0  for  transport,  2.0  for  chemistry).  The  same  input 
parameters  as  for  flame  B  are  important;  except  that,  because  of  the  low 
concentration  of  H,  the  diffusion  of  K  is  unimportant. 
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FIGURE  1.  OH  BENCHMARK  PROFILE,  FLAME  E 
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FIGURE  ?.  H9  BENCHMARK  PROFJLE,  FLAME  B 
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FIGURE  3.  NORMALIZED  LINEAR  SENSITIVITY  COEFFICIENTS  OF  THE  OH  MASS  FRACTION  WITH  RESPECT 
TO  IMPORTANT  TRANSPORT  PARAMETERS ,  FLAME  B 
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IMPORTANT  TRANSPORT  PARAMETERS,  FLAME  B 
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FIGURE  5.  NORMALIZED  LINEAR  SENSITIVITY  COEFFICIENTS  OF  THE  OH  MASS  FRACTION  WITH  RESPECT  TO 
IMPORTANT  REACTION  RATE  COEFFICIENTS,  FLAME  B 
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FIGURE  6.  NORMALIZED  LINEAR  SENSITIVITY  COEFFICIENTS  OF  THE  H0  MASS  FRACTION  WITH  RESPECT  TO 
IMPORTANT  REACTION  PATE  COEFFICIENTS,  FLAME  5 


The  sensitivity  of  the  flame  speed  to  the  input  parameters  (Tables  3  and 
4)  is  a  rough  guide  to  the  sensitivity  of  the  profiles.  The  plotted  results 
show  that  a  species  profile  is  generally  sensitive  to  a  number  of  input 
parameters.  A  similar  correlation  was  obtained  for  the  temperature  profiles 
(not  shown). 


VI.  SECOND  ORDER  SENSITIVITY  COEFFICIENTS 


In  this  section  we  develop  a  procedure  for  obtaining  second  order 
sensitivity  coefficients.  The  computations  are  rather  complicated  and  so  it 
would  neither  be  practical  nor  desirable  to  carry  them  out  for  every 
problem.  The  purpose  of  these  second  order  computations  is  to  establish 
guidelines  for  the  use  of  the  two  linear  approximations,  Equations  (18)  and 
(23).  We  do  this  by  examining  the  accuracy  and  the  region  of  validity  for 
each  of  the  linear  approximations. 

Since  the  largest  uncertainties  are  usually  associated  with  the  kinetic 
parameters,  we  limit  this  discussion  to  these  parameters.  Considering  the 
first  and  second  order  terms  of  a  Taylor  expansion  for  the  flame  speed,  we 
obtain  the  expression 
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S„  =  S  +  Z  ,  SJ(a.-l)  +  0.5  Z,  SJJ(a  -1)  + 

N  D  J-l  J  J-l  J 
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5-1  i-j+i  sJ  <v»  (“k°’ 

where 

SJk  =  Skj  =  3^S/3a.  3a,  evaluated  at  a.=  a  =  1. 

J  k  J  k 

Given  NR  rates,  there  are  NR(NR+l)/2  second  order  terras. 


(30) 


Equations  for  the  partial  derivatives,  S^,  are  found  by  taking  the 
partial  derivative  of  Eq .  (10)  and  (11)  with  respect  to  a.  and  (Kj,  k<NR) 
and  by  letting  the  a's  approach  one.  The  resulting  equations  are 
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where  the  superscripts  j  and  k  indicate  partial  derivatives  with  respect  to  a 
and  ,  respectively.  Explicitly  taking  the  partial  derivative  we  have 
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The  benchmark  values  and  the  first  order  sensitivity  coefficients  are  known, 
and  the  boundary  conditions  for  this  system  of  equations  is 


and 


1  k 

u  -  0  at  ^  =  -  00  (x=-<») 

m 

3u  3  0  at  =  +  00  (x  =  -H»). 

m 
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Thus  Eq.  (31)  and  (32)  can  be  solved  for  the  mn^  and  u^^. 


Expanding  Eq.  (24)  in  terms  of  the  lna^  UP  to  second  order,  we  obtain 
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An  equivalent  form  of  Eq. 
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As  with  the  linear  coefficients,  we  can  find  an  expression  relating  the 

second  order  sensitivity  coefficients  of  the  polynomial  expansion,  Eq.  (30) 

and  the  sensitivity  exponents  of  Eq .  (35)*  The  partial  derivative  of  Eq •  (35) 

with  respect  to  a  yields 
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By  taking  the  partial  of  Eq  .  (37)  with  respect  to  and  by  letting  the 
a's  approach  1,  we  obtain 
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Similarly,  by  taking  the  partial  with  respect  to  a^,  n  t  m  , 
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Inspection  of  Equations  (38)  and  (39)  shows  that  the  second  order 
logarithmic  coefficients,  SEmn,  are  related  to  the  second  order  normalized 
derivatives,  Smn  =  9^S/9a  9a  ,  in  a  more  complicated  fashion  than  the 

corresponding  relationshipmbetween  first  order  coefficients  given  by  Eq. 

(28).  Specifically  it  is  possible  for  SEmn  to  be  large  or  small  relative  to 
Smn/SB.  In  general,  the  expansion  (30)  in  terms  of  (a.  -  1)  and  the  expansion 
(35)  in  terms  of  In  a^  will  have  different  rates  of  convergence. 

In  the  special  case  a  =  ...  =  a  ,  the  linear  logarithmic  expansion, 
Eq.  (24)  is  exact  and  Eq.  f25)  holds.  Tnus  all  higher  order  terms  of  the 
expansion  for  InS  must  sum  to  zero.  For  the  present  case  of  a  second  order 
expansion  we  have  from  Eq.  (35)  that 
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VII  .  SECOND  ORDER  ANALYSIS  NUMERICAL  RESULTS  AND  DISCUSSION 


1.  Flame  Speed • 

The  main  difficulty  with  second  order  analysis  is  the  large  number  of 
terms.  For  example  in  the  present  case  of  17  reaction  rate  coefficients  there 
exist  153  second  order  terms.  Therefore,  in  our  examples,  second  order 
coefficients  were  computed  only  for  the  reactions  that  were  found  important  by 
the  first  order  analysis.  Table  5  lists  a  subset  of  the  flame  speed 
sensitivity  coefficients  computed  for  flame  B.  Table  h  lists  those  for  flame 
D.  Both  the  logarithmic  sensitivity  coefficients,  SE^,  and  the  normalized 
polynomial  coefficients,  S^/Sg,  are  given  for  comparison. 


TABLE  5.  SELECTED  FLAME  SPEED  FIRST  AND  SECOND  ORDER  SENSITIVITY 
COEFFICIENTS  FOR  THE  RATE  COEFFICIENTS  OF  FLAME  B. 


J 

k 

SEJ 

SEJk 

sJk/sB 

1 

1 

.0744 

-.0477 

-.1166 

2 

.0356 

.0572 

4 

.0123 

.0259 

5 

.0232 

.0314 

7 

-.0227 

-.0298 

2 

2 

.2911 

-.1211 

-.3275 

4 

.0203 

.0738 

5 

-.0361 

-.0041 

7 

.0390 

.0113 

4 

4 

.1838 

-.0493 

-.1993 

5 

.0310 

.0512 

7 

-.0248 

-.0423 

5 

5 

.1099 

-.0454 

-.1432 

7 

.0493 

.0388 

7 

7 

-.0952 

-.0624 

.0419 

0.5 

z  SEJJ 

+  E  I  S  jk  =  - 

+  Lj£k  SE 

.0359 

0.5 

E  SEJJ 

1  +  E  W  |SEJt| 

-  .4573 
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TABLE  6.  SELECTED  FLAME  SPEED  FIRST  AND  SECOND  ORDER  SENSITIVITY 
COEFFICIENTS  FOR  THE  RATE  COEFFICIENTS  OF  FLAME  D. 


j 

k 

SEj 

SEjk 

sJk/sB 

1 

1 

.0576 

-.0489 

-.1032 

2 

.0692 

.1409 

4 

-.0126 

-.0432 

5 

.0489 

.0920 

7 

-.0648 

-.1084 

2 

2 

1.2427 

-.0815 

-.4999 

4 

.6547 

-.0054 

5 

-.5327 

.3984 

7 

.6015 

-.3373 

4 

4 

-.5312 

-.7089 

.1044 

5 

.4922 

.0943 

7 

-.5160 

-.1146 

5 

5 

.7492 

-.2916 

-.4795 

7 

.3797 

-.1863 

7 

7 

-.7555 

-.4701 

.8562 

0.5  Z 

se”  + 

s  jk  =  - 

E 

.0404 

0.5  Z 

|Se3J  |  +  !  j 

2  1  S  jk 

£k  1  bE 

|  =  4.5328 

33 


are 


The  values  of  the  flame  speed  sensitivity  coefficients,  Sy  , 
negative  for  all  four  flames.  This  can  be  interpreted  by  considering  a 
simplified  case.  Consider  the  single  rate  coefficient,  j,  with  SE^>0  and 
a.  >  1.  Then  the  flame  speed  will  increase.  Physically  we  know  that,  for 
tilose  flames  that  involve  a  chain  of  reactions,  as  is  the  case  here  with 
hydrogen/oxygen,  the  faster  a  reaction  becomes  the  more  sensitive  the  flame 
speed  becomes  to  other,  slower  reactions  in  the  chain.  At  some  point  further 
increases  in  the  j  th  rate  coefficient  have  little  effect  upon  the  flame 
speed.  Numerically  then  we  expect  that  the  first  order  logarithmic 
approximation  overestimates  the  flame  speed  and  that  the  second  order 
correction  §hould  be  negative.  By  the  symmetry  of  Eq.  (23),  the  argument  also 

holds  if  Sy^<0  and  the  rate  coefficient  is  decreased,  i.e.  a.  <  1. 
h  J 

Continuing  with  our  example,  we  find  that  from  Eqs.  (18)  and  (28)  the 
linear  Taylor  expansion  in  terms  of  (cx.  -  1)  can  be  written  as 


SN(Tay)  =  SB[l+SEj  (c^  -  1)].  (41) 

The  logarithmic  formula,  Eq.  (23),  is  also  expanded  in  terms  of  (a^  -  1)  and 
the  results  are 


S  J 

SN(Log)  =  SgCtj  E  =  SB[l+SEj  (“j-O  +  S^S^-l)  (<*j-l)*  /2  +  ...].  (42) 


From  Eqs.  (41)  and  (42)  we  see  that  the  two  approximations  are  identical 
whenever  S£j  =  0  or  Sgj  =  1.  This  leads  us  to  consider  three  regions  of 
interest:  1)  0  <  Sg-^<1,  2)  Sg-^  >  1,  and  3)  Sg-^  <  0. 

The  difference  between  the  approximations  given  by  Eqs.  (41)  and  (42)  is 

SN(Log)  -  SN(Tay)  =  -  1)  (c^  -  l)2/2,  (43) 

where  we  have  assumed  that  third  order  and  higher  terms  in  Eq.  (42)  are 
negligible.  For  the  case  0  <  S£J  <  1  we  see  from  Eq.  (43)  that 

Sn(Log)  <  SN(Tay).  (44) 

We  have  argued  above  that  the  first  order  logarithmic  approximation,  SN(Log), 
will  overestimate  the  flame  speed;  then,  from  Eq.  (44),  we  expect  that  the 
approximation  given  by  Eq.  (41)  for  Sjj(Tay)  will  be  even  worse. 

We  can  check  this  heuristic  argument.  Were  the  linear  approximation  in 
terms  of  lna .  in  fact  better  than  the  linear  approximation  ip. terms  of  (a.  -  1), 
then  the  secind  order  term  in  the  logarithmic  expansion,  |SgJJ|  ought  to  be 
smaller  than  the  normalized  second  order  term  |  S^/Sgl  •  Inspection  of  the 
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values  listed  in  Tables. 5  and  6  support  this  conclusion;  i.e.  that  |  |  < 

|SJJ/SB|  whenever  0  <  S£ j  <  1 . 

For  the  cases  >  1  or  <  0  we  see  from  Eq.  (43)  that  SN(Log)  > 
S^(Tay).  In  these  cases  we  expect  the  Taylor  expansion,  Eq.  (41),  to  be  more 
accurate  than  the  logarithmic  expansion,  Eq .  (42). 

The  above  inferences  for  a  single  rate  coefficient,  can  be  summarized  and 

t 


-  D, 

Eq.  (18),  will  have  a  wider  r^nge  of  validity. 

In  addition  one  can  have  E  ~  0.5  and  E  j  S  |  large  only  if  there  are 
many  of  opposite  signs.  We  conclude  from  this  that  the  larger  the  overall 

sensitivity,  the  more  complex  the  factors  affecting  the  output  functions  are 
expected  to  be. 

In  passing  we  note  that  the  magnitude  of  the  second  order  terms  indicates 
the  degree. of  coupling  between  corresponding  reactions.  For  flame  B,  the 
largest  j*k,  is  .  Reactions  5  and  7  directly  compete  for  the 

consumption  of  HO2 ,  which  has  a  small  concentration.  For  flame  D,  the 
reactions  2,4,5  and  7  are  all  closely  coupled  because  of  the  small 
concentrations  of  H  in  this  flame.  A  change  in  any  one  of  the  above  rate 
coefficients  strongly  affects  the  other  three. 

2.  Profiles. 

We  are  also  interested  in  the  behavior  of  the  profiles.  As  examples, 
results  for  the  OH  and  H2  profiles,  flame  B,  will  again  be  used.  Figures  7 
and  8  show  the  second  order  sensitivity  coefficients:  u^^,  and 

ut44  for  OH  and  H2 ,  respectively;  the  coefficients  have  been  normalized  by 
dividing  by  the  maximum  of  each  benchmark  profile,  u^B(max).  Figures  9  and  10 
show  the  corresponding  logarithmic  coefficients.  These  are  related  to  the 
coefficients  of  the  expansion  in  terms  of  (a.  -  1);  see  Eqs.  (38)  and. (39). 

In  order  to  compare  the  two  sets  of  second  o^der  coefficients,  the  u^  are 
first  multiplied  by  u  before  normalizing  by  dividing  by  u^(max).  The 
other  sensitivity  coefficients  for  flame  B  are  not  shown  because  they  are 
small  relative  to  those  shown  in  Figures  7-10. 

Since  the  normalized  second  order  logarithmic  sensitivity  coefficients 
are  smaller  than  the  first  order  coefficients  and  since  the  second  order 
coefficients  from  the  polynomial  expansion  have  about  the  same  magnitude  as 
the  first  order  coefficients,  we  conclude  that,  in  this  example,  the 
logarithmic  formula,  Eq.  (24),  is  more  accurate  than  the  linear  polynomial 
formula,  Eq.  (18). 

There  is  no  longer  any  strong  correlation  between  large  flame  speed 
sensitivity  coefficients  and  large  profile  sensitivity  coefficients. 


with  0  <  SEJ  <  1  for  mos 
| ,  will  be  relatively 
small  and  the  logarithmic  expansion  will  be  more  accurate.  On  the  other  hand 
if  E|s  ^|is  large  (e.g.  Flame  D)  then  both  expansions  will  be  accurate  only 
for  small  departures  of  the  a.  from  unity  and  the  expansion  in  terms  of  (a. 


generalized  to  more  than  one  parameter.  For  a  flame 
SF^ ,  (e.g.  flame  B),  the  overall  sensitivity,  E | S  ^ 
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FIGURE  7.  NORMALIZED  SECOND  ORDER  SENSITIVITY  COEFFICIENTS  OF  THE  OH  .MASS  FRACTION  WITH  RESPECT 
TO  IMPORTANT  REACTION  RATE  COEFFICIENTS,  FLAME  B 
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FIGURE  0.  NORMALIZED  SECOND  ORDER  LOGARITHMIC  COEFFICIENTS  OF  THE  OH  MASS  FRACTION  WITH 
RESPECT  TO  IMPORTANT  REACTION  RATE  COEFFICIENTS,  FLAME  B 
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FIGURE  10  NORMALIZED  SECOND  ORDER  LOGARITHMIC  COEFFICIENTS  OF  THE  H2  MASS  FRACTION  WITH 
RESPECT  TO  IMPORTANT  REACTION  RATE  COEFFICIENTS,  FLAME  B 


For  flames  A  and  C  (not  shown)  the  results  are  similar  to  those  for  flame 


B. 


The  second  order  sensitivity  coefficients  of  the  OH  profile  for  flame  D 
(not  shown)  exhibit  the  following  behavior.  The  second  order  polynomial 
coefficients  are  larger  than  the  first  order  coefficients  and  the  second  order 
logarithmic  coefficients  are  even  larger.  There  are  many  coupling  terms  with 
large  values,  a  situation  similar  to  the  flame  speed  coefficients  of  Table 
6.  For  the  H2  profile  (not  shown),  the  second  order  polynomial  coefficients 
are  somewhat  smaller  than  the  first  order  coefficients  and  the  logarithmic 
coefficients  smaller  still. 

In  summary  we  find  the  results  of  these  sensitivity  analyses  follow  those 
for  the  flame  speeds.  Specifically,  for  flames  A,  B  and  C  we  find  that  the 
logarithmic  formula  is  the  more  accurate.  We  also  find  this  to  be  the  case 
for  the  major  species  and  the  temperature  of  flame  D.  However,  for  the 
radical  profiles  of  flame  D,  where  the  first  order  coefficients  are  large 
(i.e.  of  order  two),  then  both  expansions  will  be  accurate  only  for  small 
departures  of  the  a.  from  unity.  The  expansion  in  terms  of  the  (a  .  -  1)  will 
have  the  wider  range?  of  validity.  ^ 

VIII.  COMPARISONS  OF  EXTRAPOLATIONS 

So  far  we  have  developed  four  possible  extrapolations  based  on 
sensitivity  coefficients.  The  expansions  in  terms  of  (ot  .  -  1)  we  shall  call 
Linear  Taylor  and  Second  Order  Taylor.  The  expansions  irl  terras  of  lna  .  we 
shall  call  Linear  Logarithmic  and  Second  Order  Logarithmic.  In  this  section 
these  extrapolations  are  compared  with  exact  numerical  solutions  for  selected 
values  of  a  . .  The  purpose  is  to  quantify  the  heuristic  arguments  made  in  the 
last  sectiod. 

Changes  of  a  factor  of  two  in  the  input  parameters  will  be  considered. 
This  degree  of  uncertainty  is  typical  for  many  rate  coefficients,  and  is  large 
enough  that  a  linear  expansion  can  be  inaccurate. 

Tables  7  and  8  give  a  few  typical  results  for  burning  velocity  for  flame 
B.  We  consider  rate  coefficient  2,  which  has  the  largest  positive  linear 
sensitivity  coefficient,  rate  coefficient  7,  which  has  the  largest  negative 
linear  sensitivity  coefficient,  combinations  of  rate  coefficients  2  and  7,  and 
combinations  of  several  other  rate  coefficients.  The  percent  error  is 
computed  by  dividing  the  difference  of  the  numerical  value  and  the 
extrapolated  value  by  the  benchmark  value  found  in  Table  2,  and  multiplying  by 
one  hundred. 
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TABLE  7.  COMPARISONS  OF  COMPUTED  AND  EXTRAPOLATED  FLAME  SPEEDS 
FLAME  B.  LINEAR  EXPANSIONS. 


j 

“j 

Numerical 

Taylor 

Logarithmic 

% 

Errors 

2 

2.0 

299.2 

325.3 

308.3 

10.4 

3.6 

2 

0.5 

200.1 

215.3 

205.9 

6.0 

2.3 

7 

2.0 

232.3 

228.0 

235.8 

1.7 

1.4 

7 

0.5 

265.4 

263.9 

269.1 

0.6 

1.5 

2 

2.0 

311.5 

337.3 

329.3 

10.2 

7.1 

7 

0.5 

2 

0.5 

179.5 

191.3 

192.8 

4.7 

5.3 

7 

2.0 

1, 2,4,5 

2.0 

388.7 

418.0 

397.9 

11.6 

3.6 

1,2, 4, 5 

0.5 

156.2 

168.9 

159.5 

5.0 

1.3 

1,2, 4, 5 

2.0 

403.0 

430.0 

425.0 

10.7 

8.7 

7 

0.5 

1,2, 4, 5 

0.5 

141.6 

144.9 

149.3 

1.3 

3.1 

7 

2.0 
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TABLE  8.  COMPARISONS  OF  COMPUTED  AND  EXTRAPOLATED  FLAME  SPEEDS. 
FLAME  B.  SECOND  ORDER  EXPANSIONS 


j 

“j 

Numerical 

Taylor 

Logarithmic 

% 

Errors 

2 

2.0 

299.2 

284.0 

299.4 

6.0 

0.1 

2 

0.5 

200.1 

205.0 

200.0 

1.9 

0.1 

7 

2.0 

232.3 

233.2 

232.3 

0.4 

0.0 

7 

0.5 

265.4 

265.2 

265.1 

0.1 

0.1 

2 

2.0 

311.5 

295.9 

309.2 

6.2 

0.9 

7 

0.5 

2 

0.5 

179.5 

184.8 

181.0 

2.1 

0.6 

7 

2.0 

1,2, 4, 5 

2.0 

388.7 

378.3 

389.3 

4.2 

0.2 

1,2, 4, 5 

0.5 

156.2 

159.0 

156.1 

1.1 

0.1 

1, 2,4,5 

2.0 

403.0 

394.3 

401.7 

3.5 

0.5 

7 

0.5 

1,2, 4, 5 

0.5 

141.6 

143.0 

141.2 

0.6 

0.2 

7 

2.0 

All  four  extrapolations  are  accurate.  In  general,  the  linear  logarithmic 
extrapolation  is  more  accurate  than  the  linear  Taylor  extrapolation.  The  only 
exceptions  are  some  cases  involving  rate  coefficient  7.  This  is  expected  from 
the  analysis  in  Section  VII.  In  all  examples,  the  second  order  logarithmic 
extrapolation  is  more  accurate  than  the  second  Taylor  formula.  Additional 
cases  were  computed  for  flame  B  and  for  flames  A  and  C  and  the  results  are 
similar  to  these  given  here. 


Flame  D,  as  usual,  is  an  exception.  Tables  9  and  10  contain  some  results 
for  the  four  rate  coefficients  with  the  largest  sensitivity  coefficients.  The 
errors  in  the  extrapolations  are  larger  than  before,  although  still  mostly 
reasonable.  The  linear  logarithmic  extrapolation  is  slightly  more  accurate 
for  changes  in  rate  coefficient  5  (SE^  =  .7492). 
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TABLE  9.  COMPARISONS  OF  COMPUTED  AND  EXTRAPOLATED  FLAME  SPEEDS 
FLAME  D.  LINEAR  EXPANSIONS. 


j 

ot  . 

J 

Numerical 

Taylor 

Logarithmic 

% 

Errors 

2 

2.0 

23.6 

26.5 

28.0 

24.5 

36.9 

2 

0.5 

4.1 

4.5 

5.0 

3.5 

7.9 

4 

2.0 

6.6 

5.6 

8.2 

9.1 

13.2 

4 

0.5 

14.9 

15.0 

17. 1 

0.8 

18.8 

5 

2.0 

18.6 

20.7 

19.9 

18.0 

11.2 

5 

0.5 

6.6 

7.4 

7.0 

6.4 

3.4 

7 

2.0 

6.3 

2.9 

7.0 

28.9 

5.9 

7 

0.5 

18.0 

16.3 

20.0 

14.1 

17.0 

2,5 

2.0 

35.2 

43.0 

114.8 

65.8 

672.7 

4,7 

0.5 

2,5 

0.5 

0.0 

-15.2 

1.2 

128.3 

10.3 

4,7 

2.0 
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TABLE  10.  COMPARISONS  OF  COMPUTED  AND  EXTRAPOLATED  FLAME  SPEEDS 
FLAME  D.  SECOND  ORDER  EXPANSIONS. 


j  °J 

Numerical 

Taylor 

Logarithmic 

% 

Errors 

2  2.0 

23.6 

23.6 

23.1 

0.5 

4.6 

2  0.5 

4.1 

3.7 

4.1 

2.8 

0.5 

4  2.0 

6.6 

6.2 

6.9 

3.9 

2.4 

4  0.5 

14.9 

15.1 

14.4 

2.1 

3.9 

5  2.0 

18.6 

17.9 

18.5 

6.0 

0.2 

5  0.5 

6.6 

6.7 

6.6 

0.4 

0.7 

7  2.0 

6.3 

8.0 

6.3 

13.9 

0.5 

7  0.5 

18.0 

17.6 

17.8 

3.4 

1.0 

2,5  2.0 

4,7  0.5 

35.2 

45.6 

14.5 

87.5 

175.5 

2,5  0.5 

4,7  2.0 

0.0 

-8.6 

0.2 

72.2 

1.3 

For  rate  coefficient  2,  (Sg  >  1)  and  rate  coefficients  4  and  7  (Sg  <  0, 
Sg'  <  0) ,  the  linear  Taylor  extrapolation  is  usually  more  accurate  as  expected 
from  the  analysis  in  Section  VIII.  The  exception  is *  *  2.  This  change  is 

large  enough  so  that  the  asymDtotic  behavior  of  the  formulas  becomes 
important.  If  Sg^  <  0,  as  a  .  ->•  00  the  linear  Taylor  extrapolation  approaches 

•  00  while  the  logarithmic  extrapolation  approaches  0.  So  for  large  decreases 
in  the  flame  speed  (or  mass  fraction  or  temperature)  the  logarithmic 
extrapolation  will  be  more  accurate. 

When  the  above  rate  coefficients  2,5,4  and  7  are  changed,  both 
extrapolations  have  problems.  If  the  changes  are  made  so  as  to  increase  the 
flame  speed,  the  linear  Taylor  formula  is  still  reasonably  accurate.  However, 
the  second  order  approximation  is  less  accurate  than  the  first.  The  linear 
logarithmic  extrapolation  is  much  too  large.  Since  it  is  an  exponential 
function,  it  is  possible  to  badly  overestimate  for  a  set  of  large  Sg~*  s. 
However,  the  second  order  expansion  does  represent  improved  accuracy. 

If  the  changes  in  the  rates  are  made  so  as  to  decrease  the  flame  speed, 
the  asymptotic  properties  of  the  formulas  become  important.  The  numerical 
solution  indicates  extinction  of  the  flame.  The  linear  Taylor  expansion  terms 
predict  a  negative  flame  speed,  while  the  logarithmic  formula  can  produce  only 
non-negative  speeds.  So  while  the  Taylor  formula  may  be  more  accurate  for 
small  decreases  in  the  burning  velocity,  the  logarithmic  formula  is  more 
accurate  for  large  decreases. 
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Again,  we  use  flame  B  as  an  example  in  considering  the  profiles.  The 

case  a,  =  a  =  a .  =  a r  =  2.0,  a_  =  0.5,  will  be  considered. 

1  2  4  5  7 

Figure  11  shows  the  linear  extrapolations.  The  Taylor  expansion  predicts 
negative  values  near  the  cold  boundary,  which  are  cut  off  in  the  graph.  The 
logarithmic  expansion  follows  the  numerical  solution  better  here. 

At  the  first  peak  the  OH  mass  fraction  increases  substantially.  Here  the 
Taylor  series  is  more  accurate,  although  both  series  overestimate  the  correct 
value.  This  can  happen  when  a  concentration  increases  by  a  large  amount.  In 
the  flame  front  and  post  flame  region,  the  logarithmic  formula  is  more 
accurate. 

Figure  12  we  see  that  the  second  order  logarithmic  extrapolation  is  much 
more  accurate  than  the  second  order  Taylor  expansion.  In  fact,  the  Taylor 
expansion  overcorrects  badly  at  several  locations.  For  instance  near  x  =  0. 1 , 
the  negative  values  of  the  linear  extrapolations  have  become  a  spurious  OH 
peak. 


Figure  13  and  14  show  the  H2  results  over  all  the  logarithmic  expansions 
are  more  accurate. 

A  number  of  additional  cases  (not  shown)  were  run  for  flames  A,B,  and 
C.  For  the  major  species  and  temperature,  the  logarithmic  extrapolation  was 
more  accurate  for  all  the  cases  considered.  For  the  radical  species,  the 
linear  logarithmic  formula  was  usually  more  accurate  than  the  linear  terms  of 
the  Taylor  expansion,  but  there  were  exceptions.  The  second  order  logarithmic 
approximation  always  was  more  accurate  than  the  second  order  Taylor 
approximation. 

For  flame  D,  the  first  order  logarithmic  formula  for  the  radical  species 
was  very  inaccurate.  Since  radical  concentrations  were  small,  any  change  in 
the  rate  coefficients  tended  to  cause  a  large  change  in  the  radical  profiles 
(i.e.  u^jj,  >  1  ).  Even  the  second  order  logarithmic  formula  was  less  accurate 
than  the  second  order  Taylor  approximation  . 

For  the  major  species  and  temperature  profiles,  the  first  order 
logarithmic  formula  was  substantially  more  accurate  than  the  first  order 
Taylor  formula. 


XI.  SUMMARY  AMD  CONCLUSIONS 

By  employing  expansions  in  terms  of  (  a .  -  1)  and  lna .  ,  we  have 
developed  methods  and  the  corresponding  computer  programs  io  obtain  first  and 
second  order  sensitivity  coefficients.  The  codes  are  straightforward  to 
execute,  but,  for  large  sets  of  variable  parameters,  can  be  computationally 
expensive. 

We  have  applied  these  sensitivity  analyses  to  a  test  case  of  a  set  of 
H2/O2/N2  flames  that  span  a  wide  range  of  flame  speeds  and  fuel-oxidizer 
ratios.  We  found  that  the  linear  sensitivity  coefficients  for  the  flame 
speeds  by  themselves  give  a  good  overall  description  of  the  important  input 
parameters.  Conversely,  they  would  also  permit  a  given  complex  kinetic 
mechanism  to  be  understood.  (This  type  of  sensitivity  analysis  bridges  the 
gap  between  a  fixed,  determined  set  of  numbers  obtained  from  a  numerical  code 
and  the  functional  relationships  that  show  trends,  obtained  from  analytic 
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FIGURE  13.  H?  MASS  FRACTION  PROFILE,  FLAME  B:  BENCHMARK  (B) ,  NUMERICAL  SOLUTION  FOR  -  a, 

a  =  a  =  2.0,  a,  =  0.5  (N),  LINEAR  TAYLOR  EXTRAPOLATION  (T),  LINEAR  LOGARITHMIC 
4  5  / 

EXTRAPOLATION  (L) 


0.07 


II 


OJ  DC 
<3  LU 
Q 
II  DC 
O 


a 


DC 

O 


o 

o 


CO 


CO 

■ 

o 


SNOIlDVyj  ssvw 


o 

z 

CO 

o 

1 — 1 

I— 

c 

<C 

o 

-J 

i — i 

O 

DC 

CU 

LU 

C 

2: 

DC 

ZD 

1 — 

z 

X 

LU 

r\ 

DC 

CO 

O 

N _ - 

—J 

>- 

<C 

DC* 

f— 

< 

«*- 

DC 

IXI 

LU 

O 

O 

z: 

DC 

LU 

O 

CO 

O 

zc. 

•  • 

o 

CCi 

o 

LU 

LlI 

CO 

y 

C 

, — ^ 

^ - V 

u 

LU 

z: 

V — 

N _ - 

r> 

zc 

LU 

o 

_J 

LO 

1 — 1 

t — i 

• 

h- 

Lu 

o 

< 

O 

DC 

11 

o 

D- 

CU 

r-. 

cC 

>' 

<3 

DC 

O 

|— 

i— i 

9\ 

X 

1— 

O 

LU 

o 

• 

c 

C\J 

O 

DC 

►—i 

Lu 

II 

2- 

rn 

CO 

LT) 

f— 

CO 

a 

»— i 

< 

DC 

n 

c 

CO 

C\J  *3* 

o 

□c 

<3 

LU 

DC. 

ZD 

o 


49 


solutions.)  A  more  detailed  analysis  of  different  parts  of  the  flame  can  be 
had  by  studying  the  sensitivity  coefficient  profiles. 

The  adequacy  of  the  first  order  formulas  was  checked  qualitatively  by 
computing  second  order  coeffients  and  quantitatively  by  comparing  the 
expansion  predictions  with  numerical  solutions.  For  flames  with  low  overall 
sensitivity,  E  ,  the  first  order  formulas  proved  adequate  for  changes  up 

to  a  factor  of  two  in  the  variable  input  parameters.  As  the  overall 
sensitivity  increases,  the  region  of  accuracy  for  the  linear  expansions 
decreases.  We  also  found  that  if  E  S  ^  is  not  much  bigger  than  0.5  the 
logarithmic  first  order  expansion  was  the  more  accurate.  In  any  case,  the 
logarithmic  expansion  was  more  accurate  for  the  major  species  and  temperature 
profiles . 

Finally,  we  found  that  a  large  value  for  E  indicated  a  complex 

kinetic  mechanism  while  a  value  close  to  0.5  indicated  a  simple  mechanism. 
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APPENDIX  A.  INPUT  PARAMETERS 


The  enthalpies  and  heat  capacities  are  evaluated  using  the  sixth  order 
polynomial  fits  of  Gordon  and  McBride^ •  These  reproduce  the  JANAF  values 
to  within  a  few  parts  per  thousand  over  the  temperature  range  300K-3000K.  The 
constant  Cp  is  derived  from  these  values  (see  Reference  1  of  text). 

The  kinetic  scheme  consists  of  seventeen  reactions  and  has  been  taken 
from  Table  9  of  Dixon-Lewis  (see  Reference  12  of  text).  The  rates  of  the  back 
reactions  are  computed  using  the  equilibrium  constants,  which  are  calculated 
from  the  Gordon  and  McBride  polynomial  fits.  A  least  squares  fit  is  then  made 
over  the  range  300K-3000K  to  obtain  the  form  k^  =  ATB  exp  (-C/T). 

The  transport  parameters  are  given  in  Table  A-2.  The  values  for 
a  and  e /k  are  obtained  from  viscosity  measurements. 

The  viscosities  of  the  individual  species  and  binary  diffusion 
coefficients  are  calculated  as  in  Reference  1  of  the  text.  The  thermal 
conductivities  of  the  individual  species  are  computed  as  the  average  of  the 
Eucken  and  modified  Eucken  correction^  •  The  values  computed  by  this 
procedure  are  within  a  few  percent  of  the  values  computed  using  the  Mason  and 
Monchick  formula  (as  was  done  in  Reference  1  of  the  text).  This  simple 
procedure  is  used  in  this  paper  because  it  eliminates  the  necessity  of  finding 
the  rotational  collision  numbers  of  the  species,  which  are  not  well  known. 

2 

The  constants  p  D  and  pX  are  computed  from  the  above  transport 
coefficients,  using  Method  V  (see  Reference  1  of  text). 


A1S.  Gordon  and  B.J.  McBride ,  "Computer  Program  for  Calculation  of  Complex 
Chemical  Equilibrium  Compositions,  Rocket  Performance,  Incident  and 
Reflected  Shocks  and  Chapman-Jouguet  Detonations,  NASA-SP-273,  1971,  1976 
program  version . 

A^D.R.  Stull  and  H.  Prophet,  JANAF  Thermochemical  Tables,  2nd  Edition,  NSRDS- 
NBS-37,  June  1971. 

A3R.C  Reid  and  J.K.  Sherwood,  The  Praparties  of  Gases  and  Liquids,  2nd 
Edition,  McGrow  Hill,  NY.,  1966. 
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TABLE  A-l.  REACTIONS  AND  REACTION  RATE  PARAMETERS 
FOR  THE  H2-02-N2  SYSTEM 


NO 

REACTION 

A* 

B 

C 

1 

oh+h2  -  h2o+h 

1.17E09** 

1.30 

-1825 

2 

H+02  -  OH+O 

1.42E14 

0.00 

-8250 

3 

OPH2  -  OH+H 

1.80E10 

1.00 

-4480 

4 

h+o2+m'  -  ho2+m' 

1.03E18 

-0.72 

0 

5 

H+H02  -  OH+OH 

1.40E14 

0.00 

-540 

b 

f+ho2  ~  o+h2o 

1.00E13 

0.00 

-540 

7 

h+ho2  -  h2+o2 

1.25E12 

0.00 

0 

8 

oh+ho2  -  h2ok>2 

7.50E12 

0.00 

0 

9 

o+ho2  -  oh+o2 

1.40E13 

0.00 

-540 

10 

o+ho2  -  oh+o2 

1.25E12 

0.00 

0 

11 

h+h+h2  -  h2+h2 

9.20E16 

-0.60 

0 

12 

h+h+n2  -  h2+n2 

1.00E18 

-1.00 

0 

13 

w-h+o2  -  h2+o2 

1.00E18 

-1.00 

0 

14 

h+h+h2o  ~  h2+h2o 

6.00E19 

-1.25 

0 

15 

hpoh+m"  -  H20+M" 

1.60E22 

-2.00 

0 

16 

H+O+M"  -  OH+-M” 

6.20E16 

-0.60 

0 

17 

OH+OH  -0+H20 

5.7E12 

0.00 

-390 

1  = 

[h2] 

+  0.35 

[02]  +  6.5  [H20]  + 

0.44  [N2] 

1  = 

[h2] 

CM 

o 

+ 

+  5.0  [H20]  +  [N2] 

*  A  has  units  of  cm3/mole-sec  for  two  body  reactions  or  cm6/mole2-sec  for 
three  body  reactions.  The  rate  coefficient  k  is  defined  by  ATBexp  (C/T). 

**Read  1.17E09  as  1.17  x  109. 
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TABLE  A-2.  MOLECULAR  PARAMETERS  USED  FOR  THE 
DETERMINATION  OF  TRANSPORT  PROPERTIES 


SPECIES 

o[%\ 

e/k[K] 

y  [Debye] 

a[&3  ] 

References 

H 

2.050 

145.0 

0 

0 

A4 

0 

2.947 

127.2 

0 

0 

A5 

h2 

2.920 

38.0 

0 

0.79 

A4 

°2 

3.372 

128.7 

0 

1.60 

A5 

A3 

h2o 

2.600 

572.0 

1.844 

0 

A4 

n2 

3.620 

97.5 

0 

1.76 

A3 

OH  like  0 
HO2  like  O2 


A4j.  Warnatz,  "Calculation  of  the  Structure  of  Laminar  Flat  Flames  II:  Flame 
Velocity  and  Structure  of  Freely  Propagating  Hydrogen-Oxygen  and  Hydrogen-Air 
Flames ,  Ber^.  Bursan^ef  Phys_.  Chem. ,  Vol.  82,  pp.  643-649,  1978. 

A5J.M.  Heimerl  and  T.P.  Coffee,  "The  Detailed  Modeling  of  Premixed,  Laminar 
Steady-State  Flames  I.  Ozone,'  Combustion  and  Flame,  Vol.  39,  pp.  301-315, 
1980. 
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GLOSSARY 


Dim 


"fr 


br 


in 


*1 

N 

NR 

P 

R 


S 

SB 

SN 

t 

T 


_o 

concentration  of  species  i,  pY^/M^,  mole-cm 
specific  heat  of  the  mixture,  cal-gm  ^-K 

2  -1 

diffusion  coefficient  for  species  i  into  a  mixture,  cm  -s 
specific  enthalpy  for  species  i  at  temperature  TQ,  cal-gm 
specific  enthalpy  for  species  i,  equals  h°  +  Cp  (T-Tq),  cal-gm 
forward  rate  for  reaction  r,  crn  -  mole~  -s  or  cm  -mole  -s 

o  _ 1—1  f,  _ 9—1 

back  rate  for  reaction  r,  cm  -mole  -s  or  cm  -mole  -s 

mass  flux  through  the  origin  in  the  (  ij>,t)  coordinate  system, 
gm-cm~^-s  1 . 

molecular  weight  of  species  i,  gm-mole 
number  of  species, 
number  of  reactions . 
pressure,  atm. 

gas  constant  =  1.9872  cal-mole"^-K-^ . 

_3 

rate  of  production  of  species  i  by  chemical  reactions,  mole-cm 

-3  -1 

rate  of  production  of  species  i  by  reaction  r,  mole  -cm  -s 
burning  velocity  of  the  flame,  cra-s”^. 
benchmark  burning  velocity,  cm-s  ^ 

burning  velocity  for  new  values  of  the  input  parameters,  cm-s  • 
time,  s. 

Temperature,  K. 

reference  temperature  for  computing  enthalpy,  K. 
temperature  of  the  unburned  mixture,  K. 
temperature  of  the  burned  mixture,  K. 

Yi  if  i=l,2...N;  T  if  i=N+l. 
fluid  velocity,  cm  -s”^. 
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GLOSSARY 


vi>r,H  number  of  product  molecules  of  species  i  of  reaction  r 
v±9r'  number  of  reactant  molecules  of  species  i  of  reaction  r. 
x  spatial  coordinate,  cm, 

mass  fraction  of  species  i. 

Yiu  mass  fraction  of  species  i  in  the  unburned  mixture. 

Y^£  mass  fraction  of  species  i  in  the  burned  mixture, 

transformed  space  coordinate,  gm-cm^. 

6  Kronecker  delta;  equals  1  if  i-j,  equals  0  if  i  *  j. 

3-  J 

X  thermal  conductivity  of  the  mixture,  cal-cm~*-s”*-K~* . 

_o 

P  total  gas  density,  gm-cra  . 
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savings  as  far  as  man-hours/contract  dollars  saved,  operating  costs 
avoided,  efficiencies  achieved,  etc.?  If  so,  please  elaborate. 
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